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The LHCb Collaboration has recently reported the observation for the first time of a spin-3 
resonance in the heavy quark sector. They have shown that the D°K~ structure seen in the 
D°K- reaction and with invariant mass 2.86 GeV is an admixture of a spin-1 and a spin-3 
resonances. Motivated by the good agreement between our theoretical predictions some time ago and 
the properties extracted from the experiment of the D*i(2860) and D* 3 ( 2860 ) states, we perform an 
extension of the study of the strong decay properties of the D*j(2860) and present the same analysis 
for the D*i(2700) and Dsj(3040) mesons. This provides a unified and simultaneous description 
of the three higher excited charmed-strange resonances observed until now. For completeness, we 
present theoretical results for masses and strong decays of the low-lying charmed-strange mesons and 
those experimental missing states which belong to the spin-multiplets of the discovered D*i(2700), 
Db(2860) and Dsj(3040) resonances. The theoretical framework used is a constituent quark model 
which successfully describes hadron phenomenology from light to heavy quark sectors. 
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I. INTRODUCTION 

The spectrum of charmed-strange mesons contains a 
number of well established states [1] corresponding to 
the S-wave Dg and D* mesons with spin-parity 0“ and 
1“, respectively; and the P-wave states with quantum 
numbers = 0+ (Pi*o(2317)), 1+ (£>^1(2460) and 
Pisi(2536)) and 2+ (£*2(2573)). 

In addition, between the years 2006 and 2009, three 
new cs mesons were observed at the P-factories in DK 
and D*K decay modes and in three-body 6-hadron 
decays [2-4]. These states have been recently included 
in the Particle Data Group (PDG) as the £*3(2700), the 
£^^(2860) and the £sj(3040). While the £^1(2700) is 
commonly believed to have quantum numbers = 1~, 
there are several possibilities for the £*j(2860) and 
£sj(3040) states. Different predictions of the theoretical 
models can be found in Refs. [5-8]. 

Recent experiments performed by the LHGb Gollabo- 
ration have contributed to clarify the puzzle around the 
£*j(2860) resonance [9, 10]. A careful re-examination of 
the D'^K~ invariant mass around 2.86 GeV in the decay 
Bg —>■ finds that a spin-1 state and a spin-3 

state overlap under the peak. Since the resonance sub¬ 
structure of the three-body decay is analysed through a 
Dalitz plot, the well-defined initial state allows to unam¬ 
biguously determine the quantum numbers, in particular, 
the parity of the £*j(2860) state to be odd. The deter¬ 
mined masses and widths of the £*3 (2860) and £*3 (2860) 
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are [9, 10] 

M(£*3 (2860)) = (2859 ± 12 ± 6 ± 23) MeV, 
r(£*3 (2860)) = (159 ± 23 ± 27 ± 72) MeV, ^ ’ 

and 

M(£*3(2860)) = (2860.5 ± 2.6 ± 2.5 ± 6.0) MeV, 
r(£*3(2860)) = (53 ± 7 ± 4 ± 6) MeV, ^ ’ 

where the first uncertainty is statistical, the second is 
systematic and the third is due to model variations. 

Spin-3 states had never been observed in heavy 
flavoured mesons and so these new measurements have 
led to many theoretical works [11-16]. In Ref. [17] we pre¬ 
dicted two resonances 1^£3 and 1^£3 at 2.86 GeV with 
total decay widths 153 MeV and 85 MeV, respectively^. 

Motivated by the excellent agreement with the LHGb 
results, we will address in this work an extension of our 
study to all strong decay properties of the £*j(2860) 
resonance. Moreover, we will present the same analysis 
for the £*3(2700) and £5j (3040) mesons. We will 
work within the framework of a constituent quark model 
(CQM) proposed in Ref. [18] (see references [19] and [20] 
for reviews). This model successfully describes hadron 
phenomenology and hadronic reactions [21-23] and has 
recently been applied to mesons containing heavy quarks 
(see, for instance, Refs. [24-27]). 

For completeness, we will resume our theoretical 
results for the masses of the low-lying charmed-strange 
mesons and calculate strong decays for those states above 


^ Note that we are using here spectroscopic notation: n^'®+^Lj, 
where n refers to the radial excitation with n = 1 indicating the 
ground state, and 5, L and J is the spin-, angular- and total- 
momentum of the cs pair. 
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the open-flavour threshold, the -Dsi(2536) and _D*2(2573) 
mesons. There are cs states which have not yet been 
seen by experiments but belong to the spin-multiplets 
of the discovered D*;^(2700), i7*j (2860) and i7sj(3040) 
resonances. We will also compute their masses and strong 
decay properties in order to guide experimentalists in 
their search. 


This manuscript is arranged as follows. In Sec. II we 
will describe the main properties of the constituent quark 
model relevant to the heavy quark sector and review the 
^Pq strong decay model adapted to our formalism. In 
Sec. Ill we will present our theoretical results. First, we 
will review our quark model results for the properties 
of the low-lying charmed-strange mesons; second, we 
will compare the available experimental data of the 
11*^(2700), D*j{2860) and (3040) resonances with 
the theoretical predictions attending to our quantum 
number assignments. And third, we will compute masses 
and strong decay properties of those states that are 
experimentally missing and lie in the same energy range 
of the 11*^(2700), (2860) and D*j{30A0) resonances. 

Partial decay widths into all open-decay channels will 
be provided in order to guide experimentalists in the 
quest of completing the information of the charmed- 
strange meson sector. We will summarize and give some 
conclusions in Sec. IV. 


II. THEORETICAL FRAMEWORK 
A. The constituent quark model 

Spontaneous chiral symmetry breaking of the QCD 
Lagrangian together with the perturbative one-gluon 
exchange (OGE) and the nonperturbative confining 
interaction are the main pieces of constituent quark 
models. Using this idea, Vijande et al. [18] developed 
a model of the quark-quark interaction which is able 
to describe meson phenomenology from the light to the 
heavy quark sector. 

The wide energy range needed to provide a consistent 
description of light, strange and heavy mesons requires an 
effective scale-dependent strong coupling constant. We 
use the frozen coupling constant [18] 



in which ^ is the reduced mass of the qq pair and oq, 
/io and Aq are parameters of the model determined by a 
global fit to the meson spectra. 

In the heavy quark sector chiral symmetry is explicitly 
broken and Goldstone-boson exchanges do not appear. 
Thus, OGE and confinement are the only interactions 
remaining. The one-gluon exchange potential contains 
central, tensor and spin-orbit contributions given by 
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where r'o(/i) = and rg{p) = fg^^ are regulators 

which depend on the reduced mass of the qq pair. 

The contact term of the central potential has been 
regularized as 
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One characteristic of the model is the use of a screened 
linear confinement potential. This has been able to 
reproduce the degeneracy pattern observed for the higher 
excited states of light mesons [28]. As we assume that 
confining interaction is flavour independent, we hope 
that this form of the potential will be useful in our case 


because we are focusing on the high energy region of the 
charmed-strange meson spectrum. 

The different pieces of the confinement potential are 
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where Ug controls the mixture between the scalar and 
vector Lorentz structures of the confinement. At short 
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Quark masses 

m-n (MeV) 

313 


m, (MeV) 

555 


me (MeV) 

1763 


mb (MeV) 

5110 

OGE 

fo (fm) 

0.181 


fg (fm) 

0.259 


oto 

2.118 


Ao (fm“^) 

0.113 


po (MeV) 

36.976 

Confinement 

ae (MeV) 

507.4 


lie (fm“^) 

0.576 


A (MeV) 

184.432 


as 

0.81 


TABLE I. Quark model parameters. 


distances this potential presents a linear behaviour with 
an effective confinement strength a = —acfJ-d^i ' 
while it becomes constant at large distances. This type 
of potential shows a threshold defined by 

f"thr = {-ae + A}(A^Ap. (7) 

No qq bound states can be found for energies higher than 
this threshold. 

Table I shows the model parameters used herein. 
Further details about the quark model and the fine-tuned 
model parameters can be found in Refs. [18, 24, 28]. 

Among the different methods to solve the Schrddinger 
equation in order to find the quark-antiquark bound 
states, we use the Gaussian Expansion Method [29] which 
provides enough accuracy and it simplifies the subsequent 
evaluation of the decay amplitude matrix elements. 

This procedure provides the radial wave function 
solution of the Schrddinger equation as an expansion in 
terms of basis functions 


'^max 

Rdr) = C4'ni{r), ( 8 ) 

n=l I 


where a refers to the channel quantum numbers. The 
coefficients, c“, and the eigenvalue, i?, are determined 
from the Rayleigh-Ritz variational principle 

'^max 

^ (r“,„ - ENX) c“ + ^ 

n—1 L Oi' 

where and are the matrix elements of 

the kinetic energy, the normalization and the potential, 
respectively. T“,„ and fV“,„ are diagonal, whereas the 
mixing between different channels is given by 1^“ . 

Following Ref. [29] , we employ Gaussian trial functions 
with ranges in geometric progression. This enables the 
optimization of ranges employing a small number of free 
parameters. Moreover, the geometric progression is dense 
at short distances, so that it enables the description of 
the dynamics mediated by short range potentials. The 
fast damping of the Gaussian tail does not represent an 
issue, since we can choose the maximal range much longer 
than the hadronic size. 

The model described above is not able to reproduce 
the spectrum of the P-wave charmed-strange mesons. 
The inconsistency with experiment is mainly due to the 
fact that the mass splittings between the T>*g(2317), 
Z)si(2460) and iAsi(2536) mesons are not well repro¬ 
duced. The same problem appears in Lattice QGD cal¬ 
culations [30] or other quark models [31]. 

In order to improve these mass splittings we follow the 
proposal of Ref. [32] and include one-loop corrections to 
the OGE potential as derived by Gupta et al. [33]. This 
corrections shows a spin-dependent term which affects 
only mesons with different flavour quarks. 

The net result is a quark-antiquark interaction that 
can be written as: 

^dij) — diGEdij) + VcONdij) + ^OGE (1^) 

where Voge and Vcon were defined before and are 
treated non-perturbatively. Ioge°^ is the one-loop cor¬ 
rection to OGE potential which is treated perturbatively. 
As in the case of Vqge and Vcon, V"oge°*^ contains cen¬ 
tral, tensor and spin-orbit contributions given by [32] 
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where Cf = 4/3, Ca = 3, 6 o = 9, 7 b = 0.5772 and the are given in Ref. [38], here we proceed to explain briefly 


scale ^ ~ 1 GeV. 


the main ingredients in which the model is based. 


B. The ®Po Decay model 


1. Transition operator 


Meson strong decay is a complex nonperturbative 
process that has not yet been described from first 
principles of QCD. Several phenomenological models 
have been developed to deal with this topic. The most 
popular is the ^Pq model [34-36] which assumes that a 
quark-antiquark pair is created with vacuum quantum 
numbers, 

An important characteristic, apart from its simplicity, 
is that the model provides the gross features of various 
transitions with only one parameter, the strength 7 of 
the decay interaction. Some attempts have been done 
to find possible dependences of the vertex parameter 7 , 
see [37] and references therein. In Ref. [17] we performed 
a global fit to the decay widths of the mesons which 
belong to charmed, charmed-strange, hidden charm and 
hidden bottom sectors and elucidated the dependence 
on the mass scale of the ^Pq free parameter 7 . Further 
details about the global fit can be found in Ref. [17]. The 
running of the strength 7 of the ^Pq decay model is given 
by 



where fj, is the reduced mass of the quark-antiquark in 
the decaying meson and, 70 = 0.81 ± 0.02 and fio = 
(49.84 ± 2.58) MeV are parameters determined by the 
global fit. 

We get a quite reasonable global description of the 
total decay widths in all meson sectors, from light to 
heavy. All the wave functions for the mesons involved in 
the open-flavour strong decays are the solutions of the 
Schrodinger equation with the potential model described 
above and using the Gaussian Expansion Method [29]. 
Details of the resulting matrix elements for different cases 


The interaction Hamiltonian involving Dirac quark 
fields that describes the production process is given by 

77/ = a/3 5 s J (fx'ip{x)ip{x), (13) 

where we have introduced for convenience the numerical 
factor ^/3, which will be cancelled with the color factor. 

If we write the Dirac fields in second quantization and 
keep only the contribution of the interaction Hamiltonian 
which creates a (tiv) quark-antiquark pair, we arrive, af¬ 
ter a nonrelativistic reduction, to the following expression 
for the transition operator 


T = - j + Pu) 


2m, 


-•\/2®7r ; 


X 


3^1 


Pfj. -P^ 



1 al{pi,)bl{pu), 

- 0 

(14) 


where p {v) are the spin, flavour and color quantum 
numbers of the created quark (antiquark). The spin 
of the quark and antiquark is coupled to one. The 
yim{p) = P^'Yimip) is the solid harmonic defined in 
function of the spherical harmonic. 

As in Ref. [39], we fix the relation of ps with the 
dimensionless constant giving the strength of the quark- 
antiquark pair creation from the vacuum as 7 = gs/2m, 
being m the mass of the created quark (antiquark). In 
this convention, values of the scale-dependent strength 7 
in the different quark sectors following Eq. (12) can be 
found in Ref. [17]. We use herein the one corresponding 
to the charmed-strange meson sector: 7 = 0.38. 
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FIG. 1. Diagrams that can contribute to the decay width 
through the model. 


Transition amplitude 

We are interested on the transition amplitude for the 
reaction {al3)A —>■ + {^p)c- The meson A is 

formed by a quark a and antiquark /3. At some point 
it is created a {p.i') quark-antiquark pair. The created 
{pv) pair together with the (a/?) pair in the original 
meson regroups in the two outgoing mesons via a quark 
rearrangement process. These final mesons are meson 
B which is formed by the quark-antiquark pair (Je) and 


meson C with (Ap) quark-antiquark pair. 

We work in the center-of-mass reference system of 
meson A, thus we have Ka = Kq = 0 with Ka and 
Kq the total momentum of meson A and of the system 
BC with respect to a given reference system. We can 
factorize the matrix element as follow 

{BC\T\A) = S^^HKo)Ma^bc- (15) 

The initial state in second quantization is 

\A) = f d^Pad^pp5^^\KA-PA)4'A{pA)a)a[Pa)b'p{pp) |0) 

(16) 

where a (/?) are the spin, flavour and color quantum 
numbers of the quark (antiquark). The wave function 
4>a{pa) denotes a meson A in a color singlet with 
an isospin I a with projection Mj^, a total angular 
momentum Ja with projection Ma, Ja is the coupling of 
angular momentum La and spin Sa- The Pa and pp are 
the momentum of quark and antiquark, respectively. The 
Pa and pa are the total and relative momentum of the 
iajd) quark-antiquark pair within the meson A. The final 
state is more complicated than the initial one because it 
is a two-meson state. It can be written as 


\BC) = / d^Ksd^Kc {JbcMbcIitT'IJtMt) 6^^^ {K — KQ)S{k — ko) 

VI + J 

^ {JbMbJcMc\JbcMbc) {IbMij^IcMi^\IaMi^) 

Mb ,Mc,Mj^ ■'Miq 

J d^psd^p,d^pxd^ppS^^\KB - Pb)S^^HKc - Pc) 

(l)BipB)(l>c{pc)al{ps)bl{pe)alip\)blipp) | 0 ), 


(17) 


where we have followed the notation of meson A for 
the mesons B and C. We assume that the final state 
of mesons B and (7 is a spherical wave with angular 
momentum 1. The relative and total momentum of 
mesons B and C are k^ and Kq. The total spin Jbc is 
obtained coupling the total angular momentum of mesons 
B and C, and Jt is the coupling of Jbc and 1. 


3. Decay width 

The total width is the sum over the partial widths 
characterized by the quantum numbers Jbc and I 

Ta-s-bc = ^2 ^a^bc{JbcJ), (18) 

JbC:^ 


where 


The ^Pq model takes into account only diagrams 
in which the {pv) quark-antiquark pair separates into 
different final mesons. This was originally motivated 
by the experiment and it is known as the Okubo- 
Zweig-Iizuka (OZI)-rule [40-42] which tells us that the 
disconnected diagrams are more suppressed than the 
connected ones. The diagrams that can contribute to the 
decay width through the ^Pq model are shown in Fig. 1. 


rA->.Bc(TBC,/) = 27 r f dkoS{EA —EBc)\MA^Bc{ko)f ■ 

(19) 

We use relativistic phase space, so 


Ta^bc{Jbc,1) 


rriAko 


( 20 ) 
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Ds 

D* 

-‘-^s 

£>,*0(2317) 

D,i(2460) 

D,i(2536) 

D,2(2573) 

This work {as) 

1984 

2110 

2510 

2593 

2554 

2591 

This work (a^) 

1984 

2104 

2383 

2570 

2560 

2609 

Experiment 

1969.0 ± 1.4 

2112.3 ±0.5 

2318.0 ± 1.0 

2459.6 ± 0.9 

2535.18 ± 0.24 

2571.9 ±0.8 


TABLE II. Masses, in MeV, of the low-lying charmed-strange mesons predicted by the constituent quark model (os) and those 
including one-loop corrections to the OGE potential (Os). Experimental data are taken from Ref. [1]. 


where 


, _ - {ms - mc)^][m\ - {niB + mc)'^] 

^A ’ ^ ^ 

is the on-shell relative momentum of mesons B and C. 


III. RESULTS 


A. Review of low-lying states 


Table II shows the masses of the low-lying charmed- 
strange mesons predicted by the constituent quark 
model. One can see our results taking into account the 
one-gluon exchange potential (og) and including its one- 
loop corrections (Os). 

The theoretical masses for the Ds and D* mesons agree 
with the experimental measurements. We also find a 
reasonable agreement for the ZIsi(2536) and ZIs2(2573) 
masses. The state assigned to the II*q(2317) is very 
sensitive to the one-loop corrections of the OGE potential 
which bring its mass closer to the experimental one. This 
effect could explain part of its lower mass as a cs state, 
but threshold effects should be taken into account before 
a definitive statement can be given about its nature. The 
spin dependent corrections to the OGE potential are not 
enough to solve the puzzle in the I’*' cs sector. In Ref. [25] 
we have studied the = I’*' charmed-strange channel, 
finding that the ZIsi(2460) has an important non-gg 
contribution whereas the ZIsi(2536) is almost a pure gg 
state. However, the presence of non-gg degrees of freedom 
modifies the ZIsi(2536) wave function in such away that 
explains most of its decay properties [25, 43, 44]. 

Table III shows the partial and total strong decay 
widths of the mesons T>si(2536) and £>*2(2573). We 
show the absolute values in MeV and the branching 
fractions in %. One can see that the total decay widths 
reported by PDG [1] are in excellent agreement with our 
results. For the £>si(2536) meson, the PDG provides the 
following two-body decay branching ratios (concerning 
strong decays): 


Ri 


R2 


r(£)^i(2536)+ ^ D*°K+) 
r(£>^i(2536)+ ^ D*+K°) '' 

rs(£>si(2536)+ ^ D*+K°) 
r(£i«i(2536)+ ^ D*+K0) 


1.18 ±0.16, 

= 0.72 ±0.05 ±0.01, 
( 22 ) 


Meson 

n 

Channel 

(MeV) 

(%) 

r 

exp. 

(MeV) 

D,i(2536)+ 

11+ 

d*+k° 

0.43 

43.48 




D*°K+ 

0.56 

56.52 




total 

0.99 

100 

0.92 ±0.05 

D,2(2573)+ 

12+ 

D+K° 

8.02 

42.95 




D°K+ 

8.69 

46.54 




d*+k° 

0.82 

4.40 




D*°K+ 

1.06 

5.67 




£>+g 

0.08 

0.44 




total 

18.67 

100 

17 ±4 


TABLE III. Open-flavour strong decay widths, in MeV, and 
branching fractions, in %, of the Dsi(2536) and £>*2(2573) 
mesons. Experimental data are taken from Ref. [1]. 


which compare reasonably well with our theoretical 
results for i?i = 1.31 and i ?2 = 0.66. For £>*^(2573) 
meson, the PDG only reports an upper limit on the ratio 
£>*°£r+/£>°/T-l' of 0.33. Our theoretical figure 0.12 is 
compatible with such a limit. It is worth to mention 
here that in the work in which the LHCb Collaboration 
disentangles the resonance structure of the peak around 
2.86 GeV [10], they also provide the £>*2(2573) mass and 
width with significantly better precision than previous 
measurements 

M(£)^2(2573)) = (2568.39 ± 0.29 ± 0.19 ± 0.18) MeV, 
r(£i^2(2573)) = (16.9 ± 0.5 ± 0.4 ± 0.4) MeV, 

(23) 

and both are in reasonable agreement with our theoreti¬ 
cal results: 2.61 GeV and 18.67 MeV, respectively. 


B. The Dsi(2700) resonance 

It is commonly believed that the £>*]^(2700) is the first 
excitation of the £>* meson. Our quark model predicts a 
mass in this energy range (2.79 GeV) but also for the 
n'^S+iLj = 2^50 state (2.73GeV). However, if the 
£>*^(2700) had quantum numbers = 0“ it would not 
decay into DK final state, and this is incompatible with 
the experimental observations. Table IV shows the open- 
flavour strong decays of the £>*^(2700) meson as the 2^Si 
state. 
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Meson n J+ Channel 

Tapo 

B3P„ 

Texp. 


(MeV) 

(%) 

(MeV) 

DJi(2700) 21“ DK 

36.99 

21.67 


D*K 

97.78 

57.26 


DsT] 

3.67 

2.15 


Din 

9.51 

5.57 


D*K*o 

22.80 

13.35 


total 

170.75 

100 

125 ± 30 


TABLE IV. Open-flavour strong decay widths, in MeV, and 
branching fractions, in %, of the DJi(2700) meson with 
quantum numbers nJ^ =21“. Experimental data are taken 
from Ref. [1]. 


The total decay width of 11*2(2700) as the 2^S'i state 
is slightly larger but close to the experimental value. The 
information of the partial decay widths shown in Table IV 
points out that the D* K decay channel is dominant and 
the DK and D*Kq are important. In addition, the 
11 * 2 ( 2700 ) meson has traces in Dgr] and 11*77 with partial 
widths of several MeV. Finally, our theoretical value for 
the branching ratio D*K/DK is 2.6, which is a factor 3 
larger than the experimental measurement 0.91 ± 0.13 ± 
0.12 reported by the BaBar Collaboration [4]. Similar 
discrepancies can be found in other quark models. This 
fact may be an indication of a bigger mixture between the 
2^52 and l^Zl 2 states [16]. In our model the mixing is not 
fitted to the experimental data but driven by the tensor 
piece of the quark-antiquark interaction. Our states are 
almost pure ^iS '2 or ^Di in the = 1“ channel. It 
would be very helpful that the LHCb Collaboration, with 
significantly better precision than previous experiments, 
repeats the measurement of this ratio. 

C. The D*j(2860) resonance 

According to the observed decay modes, the possible 
spin-parity quantum numbers of the II*j (2860) are = 
1“, 2+, 3“, and so on. The 2+ assignment is disfavoured 
because it would be the excitation of the 11*2(2573) 
meson and our model predicts a mass around 3.1 GeV. 
Beyond J = 3 the predicted masses are much higher 
than the experimental measurement. Table V shows 
the open-flavour strong decays of the II*j(2860) as the 
third excitation of the 1“ meson (mostly dominated by 
the l^Di channel) and as the ground state of 3“ meson 
(mostly dominated by the channel). 

The predicted total decay widths are in excellent 
agreement with the LHCb observation of having two 
resonances at 2.86 GeV, one of spin-1 and another one 
of spin-3. Table V shows that the DK, D*K and DK* 
decay channels are very important for the 11 * 2 ( 2860 ) 
meson. However, the 11*3(2860) resonance decays mainly 
into DK, D*K and D*Kq final states being the partial 
width into DK* very small. We also observe that the 


Meson n 

Channel 

r^Po 


r exp. 



(MeV) 

{%) 

(MeV) 

Dlj{28m) 31“ 

DK 

53.34 

34.81 



D*K 

38.43 

25.08 



Dsn 

12.12 

7.92 



Din 

5.06 

3.30 



D*K5 

7.15 

4.67 



DK* 

37.10 

24.22 



total 

153.20 

100 

159 ± 23 ± 27 ± 72 

Dlj{28m) 13“ 

DK 

38.57 

45.32 



D*K 

26.17 

30.74 



Dsn 

1.06 

1.24 



Din 

0.35 

0.41 



D*KS 

16.16 

18.99 



DK* 

2.81 

3.30 



total 

85.12 

100 

53±7±4±6 


TABLE V. Open-flavour strong decay widths, in MeV, and 
branching fractions, in %, of the D*j(2860) meson with 
quantum numbers nj^ = 31“ or 13“. Experimental data 
are taken from Refs. [9, 10] 


11 * 2 ( 2860 ) meson has traces in DsTJ and D*t] while the 
partial decay widths of the D*^(2860) into these final 
states are very tiny. With respect the branching ratio 
measured by the BaBar Collaboration [4], we obtain 

B{D*j{2860)^D*K) jo.72 il,*2(2860), 

B{D*j{2860) ^ DK) |o.68 £>*3(2860), ^ ’ 

which compares reasonably well with the experimental 
one, 1.10 ± 0.15 ± 0.19. In view of our results, we 
cannot distinguish if the branching ratio measured by 
BaBar belongs to the 11*2(2860) or to the 11*3(2860). It 
would be again very helpful that the LHCb Collaboration 
repeats the measurement of this ratio. 

D. The Dsj(3040) resonance 

The mean 2P multiplet mass is predicted in our model 
to be 3.06 GeV which is near the mass of the Ilsj(3040) 
resonance. Therefore, the possible assignments are the 
= O’*" which only decays into DK, the I"*" which only 
decays into D*K and the 2+ which decays into DK and 
D*K. The only decay mode in which Zlsj(3040) has 
been seen until now is the D*K, and so the most possible 
assignment is that the Zlsj(3040) meson being the next 
excitation in the I"*" channel. 

Table VI shows the open-flavour strong decays of the 
Us j(3040) meson as the n = 3 1+ or 4 1+ state. The 
mass of the ilsj(3040) is large enough to allow open- 
flavour strong decays not studied before in this work 
like the D*K* final state. In fact, this decay channel is 
dominant in both n = 3 I’*' and 41+ states followed 














Meson n 

Channel 



r 

exp. 

Meson n 

Channel 

r^Pn 

&p„ 

Dsj (3040) nJ^ = 

31+ D*K 

25.22 

8.36 


D*o(2934) 2 0+ 

DK 

60.62 

32.67 


DKS 

0.76 

0.25 



Dsn 

3.31 

1.78 


D*sV 

3.26 

1.08 



D*K*o 

4.47 

2.41 


D*KS 

0.02 

0.01 



D*K* 

86.42 

46.57 


DK* 

44.28 

14.69 



DiK 

18.43 

9.93 


DtoV 

0.97 

0.32 



D[K 

0.14 

0.08 


DIK 

2.81 

0.93 



Dsn' 

12.17 

6.56 


D*K* 

156.78 

52.00 



total 

185.56 

100 


DiK 

39.81 

13.20 


75^2(3094) 2 2+ 

DK 

2.14 

0.95 


D[K 

0.69 

0.23 



D*K 

1.90 

0.84 


DIK 

11.19 

3.71 



Dsn 

0.02 

0.01 


Dscj) 

15.54 

5.15 



Din 

1.02 

0.45 


Dsi(2460)r? 

0.19 

0.07 



D*K*o 

0.80 

0.36 


total 

301.52 

100 

239 ± 35^4^ 


DK* 

40.32 

17.91 

Dsj{3040) nJ^ = 

41+ D*K 

53.48 

12.37 



D*K* 

124.26 

55.18 


DKS 

0.30 

0.07 



DiK 

17.59 

7.81 


Din 

4.97 

1.15 



D[K 

6.35 

2.82 


D*KS 

1.10 

0.25 



Dsn' 

1.71 

0.76 


DK* 

100.38 

23.21 



DIK 

26.30 

11.68 


Dion 

1.66 

0.38 



Ds(j) 

0.00 

0.00 


D*oK 

2.31 

0.53 



Dsi{2460)n 

2.44 

1.09 


D*K* 

130.91 

30.27 



D*oK*o 

0.32 

0.14 


DiK 

11.58 

2.68 



total 

225.17 

100 


D[K 

0.04 

0.01 







DIK 

123 74 

28 61 


TABLE VII. Open-flavour strong decay widths, in MeV, and 


Dscp 

1.97 

0.45 


branching fractions, in %, 
states. 

of the nJ^ = 

= 2 0+ and nj^ = 2 2+ 


Dsi{2460)n 

0.09 

0.02 







total 

432.53 

100 

239 ± 35^4^ 






TABLE VI. Open-flavour strong decay widths, in MeV, and 
branching fractions, in %, of the Dsj(3040) meson with 
quantum numbers nj^ = 31'*’or 41^. Experimental data 
are taken from Ref. [1]. 


by the DK* in the case of the 3 I'*' and by the D^K and 
DK* in the case of the 4 1+. Moreover, the n = 31+ 
state has partial decay widths in the order of tens of 
MeV for the D*K, DiK, Dgcj) and D^K decay channels; 
whereas for the nJ^ = 41+ we find only partial widths in 
the order of tens of MeV for the D*K and DiK. Finally, 
the total decay widths are large for both states, being 
that of the n J+' = 31+ state in better agreement with 
the experimental data. 

There are other two states which belong to the mean 
2P multiplet that are still missing in experiment. These 
states are the first radial excitation of the cs mesons 
with quantum numbers = 0+ and 2+, respectively. 
The masses predicted by our model are 2.93 GeV for the 
nJ^ = 2 0+ and 3.09 GeV for the 2 2+. These masses are 
in agreement with recent studies of the charmed-strange 
meson sector [13, 16]. It is worthy to remind here that the 
= 0+ channel is very sensitive to the 1-loop correction 


of the OGE potential. For the second excitation, its mass 
goes from 3.03 to 2.93 GeV. 

Table VII shows their partial and total decay widths 
calculated with the ^Pq decay model. The masses 
used for the initial mesons are the theoretical ones. 
There are less open decay channels for the nJ+’ = 
2 0+ resonance than for the 2 2+. However, the total 
decay widths are very similar being 185.55 MeV and 
225.17 MeV, respectively. The D*K* decay channel is 
dominant for the two resonances. One can also find 
traces of the DiK decay channel in both cases, but the 
DK* and D^K decays are only important for the 2 2+ 
resonance. It is remarkably that the 2 0+ resonance has 
large partial width to DK whereas this is not the case for 
the 2 2+ state. The authors of Ref. [16] predict the same 
behaviour than us and notice that there is an evidence of 
a structure around 2.96 GeV in the D^K~ invariant mass 
spectrum given by LHGb [9, 10] that can be associated 
to the 2 0+ resonance. 


E. The missing states around 2.86 GeV 

We have assigned to the 11*2(2700) meson the nJ^ = 
21 “ state and to the Il*j(2860) meson the 31“ and 
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Meson n 

Channel 

Tapn 


Dsi2729) 2 0“ 

D*K 

190.62 

66.81 


DKS 

87.41 

30.64 


D*sV 

7.28 

2.55 


total 

285.31 

100 

Ds2{2888) 12“ 

D*K 

54.74 

23.45 


DKS 

35.92 

15.39 


DSv 

0.99 

0.43 


D*KS 

8.19 

3.51 


DK* 

133.55 

57.22 


DSoV 

0.00 

0.00 


DSK 

0.00 

0.00 


total 

233.39 

100 

L>„2(2948) 2 2“ 

D*K 

85.70 

45.86 


DKS 

4.52 

2.42 


DSv 

19.10 

10.22 


D*KS 

9.51 

5.09 


DK* 

31.61 

16.92 


DSoV 

0.00 

0.00 


DSK 

0.02 

0.01 


D*K* 

36.30 

19.42 


DiK 

0.12 

0.06 


D[K 

0.01 

0.00 


DsTj' 

0.00 

0.00 


total 

186.89 

100 


TABLE VIII. Open-flavour strong decay widths, in MeV, and 
branching fractions, in %, of the nj^ =20“, nJ^ = 12“ and 
nJ^ =2 2“ states. 


13“ states. The last two assignments are based on 
the disentanglement of spin-1 and spin-3 resonances 
around the 2.86 GeV peak performed by the LHCb 
Collaboration [9, 10]. Around this energy range, there are 
still three states not seen by experiments. These are the 
second radial excitation of the charmed-strange ground 
state n = 2 0“ and the 1 2“ and 2 2“ states which are 
the H-wave partners of the D*j(2860) meson. 

The masses predicted by our constituent quark model 
are 2.79 GeV, 2.89 GeV and 2.95 GeV for the nJ^ = 
2 0“, 12“ and 2 2“, respectively. These masses are in 
agreement with Refs. [13, 16] except for the 2 0“ state 
which seems that our model predicts a mass slightly 
higher. Table VIII shows partial and total decay widths 
for the three states. One can see that we predict a 
2 0“ state much broader than recent studies of the same 
states [13, 16]. There are two main reasons for this: i) 
our theoretical mass is different and this influences the 
calculation of the decay widths in the model; ii) 
Refs. [13, 16] do not calculate the partial width into the 
FKq despite this decay channel is open for the mass they 
predict. The FKq decay channel contributes 30% to the 
total decay width in our model and sum almost 100 MeV 
to it. 


With respect the other two states, 1 2“ and 2 2“, their 
masses and total decay widths are very similar to those 
predicted in Refs. [13, 16]. One can distinguish common 
features between their predictions and ours. The first 
one is that these states seems quite broad with total 
decay widths around 200 MeV. The second one is that 
the dominant decay channels for the 2 2“ state are D*K, 
DK* and D*K* with traces also in the D*ri. However, 
there are also differences. The relative order in the 
dominant channels for the 2 2“ state is different and, 
what is more important, our predictions for the partial 
decay widths of the nJ^ = 12“ seems quite different. 
They predict that its dominant decay channel is D*K 
followed by the 0^7] whereas we have a dominant DK* 
followed by the D*K and DKq decay channels. 

IV. SUMMARY AND CONCLUSIONS 

We have performed an extensive study of strong decay 
properties for the 11*^(2860) and 11*3(2860) resonances. 
We have completed the study with the same analysis for 
the 11*3(2700) and Ilsj(3040). Our theoretical results 
indicate that the 11*3(2700), 11*3(2860), 11*3(2860) and 
Zlsj(3040) mesons can be accommodated as the nJ^ = 
2 1“, 3 1“, 1 3“ and 3 I'*', respectively. These predictions 
are in agreement with other theoretical studies of the 
same resonances. 

For this study we have used a constituent quark model 
which describes successfully the hadron phenomenology 
from light to heavy quark sectors. The 1-loop corrections 
to the OGE potential and the coupling of non-qq degrees 
of freedom in the I"*" cs channel improve the theoretical 
mass splittings between the Zl*o(2317), Zlsi(2460) and 
Zlsi(2536) mesons. Moreover, the decay properties of 
the Ils3(2536) and Z1 ,j 2(2573) are well reproduced. It 
is worth to mention that the strong decays have been 
calculated using an adapted version of the ^Pq decay 
model in which the strength 7 of the decay interaction 
depends on the mass scale through the reduced mass of 
the quark-antiquark in the decaying meson. 

Finally, there are states still undiscovered by exper¬ 
iments in the mass energy region of the 11*3(2700), 
Zl*j(2860) and Zlsj(3040) mesons. We have provided 
the masses and strong decay properties of those which 
belong to the same spin-multiplets. We hope that this 
study will help experimentalists in carrying out a search 
for them. 
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